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ABSTRACT 

The first high resolution Spitzer IRS 9-37/ini spectra of 29 Seyfert galaxies (about one quarter) of 
the 12/ini Active Galaxy Sample are presented and discussed. The high resolution spectroscopy was 
obtained with corresponding off-source observations. This allows excellent background subtraction, 
so that the continuum levels and strengths of weak emission lines are accurately measured. The result 
is several new combinations of emission line ratios, line/continuum and continuum/continuum ratios 
that turn out to be effective diagnostics of the strength of the AGN component in the IR emission of 
these galaxies. The line ratios [NeV] / [Nell] , [OIV]/[NeII], aheady known, but also [NeIII]/[NeII] and 
[NeV]/[SiII] can all be effectively used to measure the dominance of the AGN. We extend the analysis, 
already done using the 6.2/im PAH emission feature, to the equivalent width of the 11.25/im PAH 
feature, which also anti-correlates with the dominance of the AGN. We measure that the 11.25/^m 
PAH feature has a constant ratio with the H2 S(l) irrespective of Seyfert type, approximately 10 to 
1. Using the ratio of accurate flux measurements at about 19/im with the two spectrometer channels, 
having aperture areas differing by a factor 4, we measured the source extendness and correlated it 
with the emission line and PAH feature equivalent widths. The extendness of the source gives another 
measure of the AGN dominance and correlates both with the EWs of [Nell] and PAH emission. Using 
the rotational transitions of H2 we were able to estimate temperatures (200-300K) and masses (1-10 
X 10^ Mq), or significant limits on them, for the warm molecular component in the galaxies observed. 
Finally we used the ratios of the doublets of [NeV] and of [SIII] to estimate the gas electron density, 
which appears to be of the order of Ue ~ 10"^^^ cm^"^ for the highly ionized component and a factor 
10 lower for the intermediate ionization gas. 

Subject headings: Galaxies: Active - Galaxies: Starbursts - Infrared: Galaxies 



1. INTRODUCTION 

Mid-infrared (mid-IR) spectroscopy provides a pow- 
erful tool to investigate the nature and physical pro- 
cesses in active galactic nuclei (AGNs) and in the star- 
burst dominated regions frequently associated to them. 
Because of the large variety of fine structure lines 
present in the mid-IR, covering a wide range in ioniza- 
tion/excitation conditions and gas density (e.g. Spinoglio 
& Malkan 1992) mid-IR spectroscopy of the Narrow Line 
Regions (NLR) in AGNs can add information not avail- 
able from classical optical spectroscopy, especially when 
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dust extinction is high. Furthermore, the electronic tran- 
sitions responsible for the infrared emission lines of vari- 
ous elements are less sensitive to uncertainties in temper- 
ature than the corresponding optical lines. Moreover the 
brightest H2 rotational lines, that can be used to quan- 
tify the presence and excitation of warm molecular g as, as 
well as a prominent Polycyclic Aromatic Hydrocarbons 
(Puget & Leger 1989), hereafter PAH, feature at 11.25 
/xm are in the mid-IR. The first detailed mid-IR spec- 
troscopic studies of AGNs and Starburst galaxies using 
multiple ionic transitions of various elements were per- 
formed by various authors (Sturm et al. 2002; Spinoglio 
et al. 2005; Verma et al. 2003, 2005, for a review) with 
the Short Wavelength Spectrometer (SWS) (de Graauw 
et al. 1996) onboard of the Infrared Space Observatory 
(ISO) (Kessler et al. 1996). However, the improved sen- 
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sitivity of the Infrared Spectrometer (IRS) (Houck et al. 
2004) on the Spitzer Space Telescope now enables a de- 
tailed investigation of the nature and physical processes 
in large samples of galactic nuclei from nonstellar (AGN) 
and stellar (starburst) power sources. 

Even at low resolution, the IRS spectra of the first few 
classical AGNs already showed the diversity of the mid- 
IR spectral features: silicate absorption and emission, 
PAH emission and strong fine structure lines (Weedman 
et al. 2005). 

Buchanan et al. (2006) examined 51 low resolution 
IRS spectra of 12/im selected Seyfert galaxies (Rush, 
Malkan & Spinogho 1993), exactly the sample we are 
considering in our study. They report a few major find- 
ings: (1) the sample contains a very wide range of con- 
tinuum types, with no more than about 3 galaxies being 
closely similar to one another, however principal com- 
ponent analysis applied to their data suggests that the 
relative contribution of starburst emission may be the 
dominant cause of variance in the observed spectra; (2) 
the starburst component in the sample objects does not 
contribute more than 40% of the total IR flux density; (3) 
Seyfert I's have higher ratios of infrared to radio emis- 
sion (sec also Rush, Malkan & Edelson 1996); (4) the 
Seyfert 2 galaxies typically show stronger starburst con- 
tributions than Seyfert I's. 

Gorjian et al (2007) foimd a strong correlation be- 
tween the [NeV]14.3/im and the [NeIII]15.5/xm lines in 
Narrow Line Regions (NLR) of AGNs, spanning 4 orders 
of magnitude in luminosity. This would imply a very nar- 
row range in ionization parameter (-2.8<logU<-2.5) for 
simple constant density photoionization models. Dudik 
et al. (2007) discussed the ratio of the doublets of [NeV] 
and [SIII] in a heterogeneous sample of active galaxies. 
Finally ULIRG galaxies, more than Seyfert galaxies, have 
been so far the object of systematic spectroscopic stud- 
ies with Spitzer IRS (Higdon et al. 2006; Armus et al. 
2007; Desai et al. 2007; Farrah et al. 2007). 

In this article we present the first high resolution IRS 
spectra of 29 galaxies from the original list of Seyfert 
galaxies of the " 12/zm Galaxy Sample" (12MGS) (Rush, 
Malkan & Spinoglio 1993), an IRAS-selected all-sky sur- 
vey flux-limited to 0.22 Jy at 12/im. 

The sample selection is briefly described in §2, the ob- 
servations are described in §3, the results are presented 
in §4, the diagnostics diagrams using line ratios, equiva- 
lent widths and other observed quantities are presented 
and discussed in §5 and a comparison between Seyfert 
l"s and Seyfert 2's is given in §6. The conclusions are 
given in §7. 

2. THE SEYFERT GALAXIES OF THE 12/iM GALAXY 
SAMPLE 

This sample is essentially a bolometric flux-limited sur- 
vey outside the galactic plane, because of the empir- 
ical fact that all galaxies emit a constant fraction of 
their total bolometric luminosity at 12/im. This frac- 
tion is - 15% for AGNs (Spinogho & Malkan 1989) 
and ^ 7%} for normal and starburst galaxies, indepen- 
dent of star formation activity (Spinoglio et al. 1995). 
Moreover, the 12MGS is less subject to contamination 
by high star-formation rate objects than other infrared 
samples defined at longer wavelengths (Hunt & Malkan 
1999). The 12MGS contains: 53 Seyfert Is and 63 Seyfert 



2s (Rush, Malkan & Spinoglio 1993). A unique ad- 
vantage of our sample is that it already has the best 
and most complete set of observations at virtually every 
other wavelength: full IRAS and near-IR coverage (Rush, 
Malkan & Spinoglio 1993; Spinoglio et al. 1995), X rays 
(Rush et al. 1996), optical spectroscopy, radio (Rush, 
Malkan & Edelson 1996), optical/IR imaging (Hunt et 
al. 1999; Hunt & Malkan 1999), ISOPHOT data on half 
of the Seyfert galaxies (Spinoglio, Andreani & Malkan 
2002). In the recent years lO/xm imaging (Gorjian et al. 
2004), K band and 2.8-4.1/xm slit spectroscopy (Iman- 
ishi 2003; Imanishi & Alonso-Herrcro 2004) optical spec- 
tropolarimetry (Tran 2001, 2003) and radio observations 
(Thean et al 2000, 2001) have been collected for most of 
the Seyfert galaxies in our sample. Low resolution IRS 
Spitzer spectra have been collected on 87 Seyfert galaxies 
of the 12/im sample and presented for a large fraction of 
them (51 objects) (Buchanan et al. 2006). 

3. OBSERVATIONS AND DATA REDUCTION 

The Spitzer IRS observations presented here have been 
collected with the Spitzer Cycle 3 GTO program "IRS 
spectroscopy of a complete sample of Seyfert galaxies in 
the Local Universe" (P30291, PI. G. Fazio). In following 
papers we will present the remaining data of this program 
as well as the results of a complete reduction of the other 
high resolution spectra of the 12/im Seyfert galaxies of 
the Spitzer archive. We refer to another future paper for 
the comparison of the data with photoionization models. 

A complete description of the IRS onboard of Spitzer 
Space Telescope can be found in Houck et al. (2004) . For 
our observations we used the two high resolution modules 
of IRS, that allow observations at the wavelength range 
from 10 to 37/im with resolutions R=A/(5A=600. The 
short wavelengths in the range 9.7-19.6 /xm are covered 
by the Short-High (SH) module and the 18.7-37.2 /im 
range by Long- High (LH). 

We acquired the targets using the red (22/im) IRS 
peak-up camera in high accuracy mode to locate the 
mid-IR centroid of the source and then offset to the ap- 
propriate slit using the standard IRS staring observing 
mode. For each observation of a Seyfert nucleus, we also 
collected an observation of an off-source position at a 
distance of 2' in declination, to be able to subtract the 
background emission from each on-source observation. 
To allow for redundancy in the data, each source was 
observed multiple times (cycles) at two positions (the 
nods) along each slit at 1/3 and 2/3 of its length, both 
with the requested integration time. 

In Table 1 the journal of the observations is presented, 
which gives for each source its equatorial coordinates, the 
Seyfert type, the redshift, the date of the observation and 
the total on-source integration times in the SH and LH 
slits. As it can be seen from the table, for every observa- 
tion, we have obtained multiple cycles so that each of the 
two nod positions of each spectrograph (SH - LH), both 
on and off-source position, arc sampled multiple times. 
This not only increases the signal to noise ratio but is 
also essential to correct for glitches in the pixel response, 
mainly due to cosmic rays hits. 

The data were processed by the Spitzer Science Center 
(SSC) pipeline version 15. The three-dimensional data 
cubes were converted to two-dimensional slope images 
after linearization correction, subtraction of darks, and 
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cosmic-ray removal. To rcdiicc the data, we used the 
software package Spectroscopy Modeling Analysis and 
Reduction Tool (SMART) (Higdon et al. 2004), specif- 
ically designed for the reduction of IRS spectra. The 
starting files are the basic calibrated data (BCD), which 
are provided by the SSC. The first step in the data re- 
duction consisted in removing the bad-responsing pixels 
by cleaning the BCD images with the rogue pixel masks 
of the IRS campaign relative to each observation. 

Off-source images were simply averaged to obtain a 
reliable background image at each nod position of each 
spectrograph. The extraction of the spectra from the two 
dimensional cleaned BCD images was done as " full aper- 
ture extraction" (see the SMART manual) , using the off- 
source averaged spectra as "sky images" . The resulting 
spectra of each nod position of each spectrograph were 
compared to remove the border effects of each order and 
the glitches. As many cleaned spectra as cycles were ob- 
tained for each nod position. Finally, the cleaned spectra 
from the same nod were averaged, using weighted means 
with 2.5(7 clipping, to get one final single cleaned spec- 
trum. 

We note that three galaxies in Table 1 , which were clas- 
sified as Seyfert galaxies at the time of the 12/zm sample 
selection (Rush, Malkan & Spinoglio 1993), have been re- 
classified (Tran 2003) as LINERs or startburst galaxies: 
these are: MCG-03-34-63, NGC7496, NGC7590. The fi- 
nal Seyfert list of this study contains therefore 13 type 
I's (including both galaxies of the pair of MRK1034) and 
14 Type 2's. We therefore use a different symbol for the 3 
non-Seyfert galaxies in the plots and exclude them from 
any analysis comparing type I's with type 2's. 

4. RESULTS 

The only emission feature which is spectroscopically 
resolved is the 11.25/im PAH blend, while the fine struc- 
ture and H2 emission lines are not resolved with IRS. 
The line fluxes have been measured using Gaussian fits, 
done within the SMART software package. In Table 2 we 
report the line fluxes, together with their la statistical 
errors, of the common ionic fine structure lines, while we 
list in Table Al, in the Appendix, the line fluxes of weak 
lines that have been detected only on a few galaxies. In 
Table 3 we list the fluxes and Icr errors of the pure rota- 
tional transitions of II2 and the PAH integrated fluxes, 
measured with the moment fit, as well as the PAH equiv- 
alent widths. When the signal to noise ratio was less than 
3, we have quoted in the tables 3a upper limits. 

The PAH emission feature is detected in all sources 
except for IRAS00521-7054, where the red part of the 
silicate absorption feature could hide the PAH emission. 
Figures 1 show the measured high resolution IRS spectra. 

4.1. Source extendness 

Although the continua have flux densities generally 
rising with wavelength, the larger aperture for the 
LH spectra results in some cases in a substantially 
higher continuum than measured at the red end of the 
SH spectrum, especially in IRAS00198-7926, NGC424 
(with major and minor diameters: 1.8'x0.8'), NGC513 
(0.7'x0.3'), NGC931 (3.9'x0.8'), NGC1125 (1.8'x0.9'), 
NGC4501 (6.9'x3.7'), NGC4748, MKN897 (0.6'x0.6') 
andNGC7590 (2.7' x 1.0'). This is due to the difference in 
the two entrance slits of the spectrometers (4". 7x11". 3 



and 11". 1x22". 3 for SH and LH, respectively). This in- 
crease in area covered by a factor of 4 has only minor 
effects when the source is more compact. In particular, 
the nonstellar AGN continuum is probably unresolved 
by even the smaller aperture (Edelson & Malkan 1986; 
Edelson et al. 1987). 

Because our spectra have been carefully background- 
subtracted, and there is no contamination from zodiacal 
emission in our slits, by measuring the flux values at the 
longest wavelength of the SH spectrograph and at the 
shortest one of the LH, we can use the ratio of these two 
fluxes as a measure of how much each of our galaxies de- 
viates from an unresolved point source. We call the ratio 
between the larger (LH) to the smaller (SH) aperture 
fluxes extendness factor, exactly the inverse of the "com- 
pactness" ratio defined by Edelson & Malkan (1986); 
Edelson et al. (1987), which is the ratio of small beam 
to large-beam 12/xm photometry. We have reported in 
table 4 the photometric fluxes measured at 19/im in SH 
(using a bandwidth of about 1 /iin) and the fluxes mea- 
sured at 19.5/um in LH (in a similar bandwidth), exclud- 
ing the [SIII]18.71/im line and the additional noise fre- 
quently present at the border of the orders. The error 
given is the standard deviation of the linear fit we per- 
formed. We divided our sources in three classes: I) the 
compact sources with R< 1.15 (10 galaxies); II) moder- 
ately extended sources with 1.15<R<1.5 (14 galaxies); 
and HI) very extended sources with R> 1.5 (6 sources). 
We allow for compact sources an increase of 10% in flux 
from the SH and the LH waveband, because all the spec- 
tra of our targets are moderately rising with wavelength. 

4.2. H2 excitation diagrams: H2 temperatures and 

masses 

From the measured H2 rotational line fluxes, we de- 
rived in Figures 2 the excitation diagrams, using the 
Boltzmann equation, which relates the level population 
to the temperature, following the method described by 
Rigopoulou et al. (2002). In Table 5 we report the H2 
temperatures, as derived from the excitation diagrams 
of Fig. 2. We estimated the temperatures, or their up- 
per or lower limits, from the slopes obtained connecting 
two adjacent measurements in the excitation diagrams, 
namely the T(O-l) from the S(0)-S(1) fit, the T(l-2) from 
the S(l)-S(2) and the T(2-3) from the S(2)-S(3). 

From the H2 line column densities, following 
Rigopoulou et al. (2002), we then estimated the warm 
H2 masses, considering the excitation temperature as de- 
rived from the S(0)-S(1) fit and the S(l) line flux. The 
estimated mass is related to the IRS aperture and it as- 
sumes that the S(l) emission fills the beam. If the galaxy 
is more extended the masses derived have to be consid- 
ered as lower limits, on the contrary, if the H2 emission is 
more compact than the IRS SH beam, the derived num- 
bers have to be considered upper limits. 

For eight galaxies we can derive an average mean tem- 
perature from the temperature upper and lower limits 
given in Table 5: MRK1034 NED02 (T - 188K), MRK6 
(T - 246K), NGC3516 (T - 263K), NGC4501 (T - 
300K), NGC4748 and NGC4968 (T ~ 250K), MCG-03- 
34-63 ( T ~ 217K), NGC7496 (T ~ 240K). For NGC513, 
the pair of detected lines infer a temperature of 191K. 
Five galaxies show a much flatter slope between the S(2) 
to S(3) lines or S(l) to S(2) than between the S(0) to 
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S(l) lines: these are ESO545-G013, MRK9, NGC4602, 
IRASF22017+0319 and NGC7590: this reflects distinct 
temperature components in the gas: a warmer compo- 
nent of about or more than 400K and a cooler one of 
about or less than 300K. For the remaining 15 galaxies 
only a broad range in temperature, as wide as 100-150K 
or more, can be given. 

There is no obvious correlation of the excitation tem- 
peratures of H2 with Seyfert type or any other properties 
in Table 2. 

The estimated H2 masses are significantly smaller than 
those found by the same procedure in ultra-luminous in- 
frared galaxies (Higdon et al. 2006), but arc comparable 
to those observed in more normal spiral galaxies (Roussel 
et al. 2007). Thus there is no indication of an unusu- 
ally large reservoir of molecular gas in the centers of our 
Seyfert galaxies. 

5. AGN DIAGNOSTIC DIAGRAMS 

5.1. Line ratios versus PAH equivalent widths 

Ratios of emission lines of different ionization poten- 
tial for the same element have been used to determine 
the hardness of the underlying radiation field in galax- 
ies (for example see Spinoglio & Malkan 1992; Verma 
et al. 2003). This allows us to quantify how much 
IR emission arises from O and B stars, and how much 
from an AGN, which has much stronger EUV and X- 
ray emission. The mid-IR lines of neon [NeII]12.81y^m, 
[NeIII]15.56/im, [NcV]14.32/im, with ionization poten- 
tials of 21.6eV, 41.0eV and 97.1eV, are among the prime 
candidates for tracing emissing from AGN (e.g. Sturm 
et al. 2002). Their proximity in wavelength sugg ests 
that effects of differential extinction are minimized. In 
addition the high ionization potential of [NeV] ensures 
that it has been excited by photons originating from the 
accretion disk of an AGN. Since the [NeV] line is rather 
faint, in particular for the ISO SWS spectrometers the 
[OIV]25.91/xm line with an ionization potential of 54.9eV 
has also been used (e.g. Genzel et al. 1998; Armus et 
al. 2007; Farrah et al. 2007). We note though that 
[OIV]25.91/xm is not as clear a diagnostic since it could 
also be excited by shocks, or extremely metal-poor star 
clusters (Lutz 2002). 

With ISO spectra of ULIRGS, Genzel et al. (1998) 
used the inverse correlation of [NeV] /[Nell] and 
[OIV]/[NeII] with equivalent width of the 6.2/zm PAH 
emission to infer the proportion of IR emission produced 
by the central AGN. Recently, Desai et al. (2007) have 
shown that the EW of the 6.2/xm PAH emission is sys- 
tematically weaker in the most luminous ULIRGs, partic- 
ularly those with warm IRAS colors. The PAH EW can 
be small even if their optical spectra show no hint of an 
AGN. They suggest that these galaxies may host strong 
star formation as well as obscured AGN activity. Instead 
of the PAH 6.2/im emission feature, commonly used in 
low resolution shorter wavelength spectra, our data al- 
lowed an accurate measurement of the 11.25/xm PAH 
emission feature, which is 0.14 % of the Lir (Soifer et 
al. 2002). We then use the [NeV]14.3/im/[NeII]12.8^m 



line ratio, the [OIV 
tio, the [NeV]14.3/xm/' 



14.3/im/[NeII]12.8/um line ra- 
SiII]34.8/um line ratio and the 

[NeIII]15.5/xm/[NeII]12.8/im line ratio as a function of 
the PAH 11.25/tim equivalent width in the diagrams of 
Fig. 3 (a, b) and Fig. 4 (a, b), respectively. 



Throughout this paper, we have looked at the corre- 
lations using the following method: first we applied to 
the data a weighted least squares (hereafter w.l.s.) fit, 
to derive a slope for each data set, considering both the 
Seyfert I's and 2's separately and all the Seyfert's to- 
gether. To estimate the significance of the correlation 
found we have then applied a bootstrap method (with a 
1000 times resampling). From each bootstrapped sample 
we have computed the slope with a w.l.s. fit and then 
used the mean value of these slopes and the standard de- 
viation. We considered valid a correlation if the resulting 
slopes from the two methods, i.e. simple w.l.s. fit and 
bootstrapped w.l.s. fit, were consistent within a thresh- 
old of ± 0.25. For each fit where a correlation is found, 
we used the mean value of the slope, as derived from the 
bootstrap method, and its standard deviation to give a 
confidence interval. 

For all the Seyfert galaxies, but especially 
for Seyfert 2's, it appears from Fig. 3 and 
Fig. 4 that there is a clear trend between the 
[NeV]14.3Aim/[NeII]12.8Aim, [NeV]14.3Aim/[SiII]34.8Aim 
and [NeIII]15.5/im/[NeII]12.8/im line ratios and the PAH 
equivalent width. The ratio of the high-ionization lines, 
excited by the AGN, to the low-ionization line excited 
in HII regions ([NeII]12.8/im) and also to [SiII]34.8/im, 
increases by about two orders of magnitude as the PAH 
equivalent width decreases. However we do not find any 
significant correlation. 

Because of the poor statistics, we cannot exclude that 
such correlations could appear with a larger sample of 
objects, and we postpone to the next articles, which will 
present new spectra, any firm conclusion. 

5.2. Line Equivalent Width Diagnostics 

As with the ratios of high- to low-ionization lines, the 
equivalent width can also function in the same way to 
measure AGN dominance. The ratio of an H Il-region 
dominated line to the underlying continuum dominated 
by the AGN should decrease systematically as the AGN 
produces an increasing proportion of the IR luminosity. 
In Fig. 5 the [NeII]12.8/xm equivalent width as a fuction 
of the PAH equivalent width is presented. The equiva- 
lent width of the PAH shows an increasing trend with 
the [Nell] equivalent width, however no correlation is 
present. 

5.2.1. Line equivalent widths versus extendness factor 

Edelson & Malkan (1986) and Edelson et al. (1987) 
first considered the extendness of the CfA Seyfert galax- 
ies as a useful AGN measure, using the ratio of 10/xm 
ground-based small-aperture flux to the IRAS 12(Um flux. 
We have correlated the extendness factor at 19/xm de- 
fined in §4.1 and listed in Table 4 with the equivalent 
widths of the key lines measured in our spectra (listed 
in Table A2 in the Appendix) and the PAH feature. We 
present the results in figures 6-8. In particular we show 
in Fig 6 (a and b, respectively) the EW of the 11.25/im 
emission feature and that of the [NeII]12.8/^m line against 
the; cxtcindness factor. We can sec that the trend in these 
two diagrams is the same: increasing the size of the 19/zm 
emitting region, both the PAH and the [Nell] equivalent 
widths increase with a very similar relationship. This 
confirms that the amount of star formation, given by 
the integrated number of HII regions due to young stars. 
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which is measured by the [Nell] Une, or by the strength 
of the PAH emission over the continuum nicely correlates 
with the source extendness. For a compact source both 
emissions are minimized, because of the dominance of 
the point-like AGN. 

In Fig. 7 the EW of the [NeIII]15.5^m line against the 
extendness factor is presented. This diagram does not 
show a strong correlation as the previous ones, because 
the [Nelll] emission lines is not simply related to the 
amount of star formation/young stars, but it includes 
contributions from the AGN: being photoionized by the 
AGN radiation field in the NLR. 

In Fig. 8 (a and b, respectively) the [NeV]14.3/im and 
[OIV]25.9/im line equivalent widths as a function of the 
extendness parameter are presented. There is no correla- 
tion left because, as we expect, the emission in these lines 
is mostly due to the AGN, just like most of the under- 
lying continuum. With a significant scatter, the average 
equivalent widths of the [NeV] and [OIV] lines of 0.02 
and 0.07/im represent "pure" AGN emission. 

5.2.2. Line ratios as a function of luminosity and 
spectral index 

One of the reasons why the 12/Ltm band has been used 
to select a complete sample of Seyfert galaxes in the lo- 
cal universe was because this band is carrying a constant 
fraction of the bolometric luminosity of Seyfert galaxies 
(Spinoglio & Malkan 1989). We consider therefore that 
for our sample the 12/xm luminosity is a good approxi- 
mation of the total luminosity and use this to see correla- 
tions with the power of the AGN, as indicated by the line 
ratio of [NeV] /[Nell] and to a lesser extent by the ratio 
[NeIII]/[NeII]. We show the [NeV]14.3Aim/[NeII]12.8Atm 
line ratio versus 12y^m luminosity in Fig. 9. Applying 
our statistical method, no correlation is found. 

Another method that has been used with the large- 
beam wideband photometry (e.g. IRAS), followed the 
discovery that local Seyfert galaxies, and AGNs in gen- 
eral at low redshift, have a "warm" mid-to-far-IR spec- 
tral index (see, e.g. Rush, Malkan & Spinoglio 1993, for 
the 12/im selected sample). We show in Fig. 10 (a and 
b, respectively) how the 60-25/im spectral index fiattens 
as the two ionization-sensitive line ratios [NeV] /[Nell] 
and [NeIII]/[NeII] increase. The higher these line ratios, 
the flatter -or warmer- is the spectral index. Because a 
warmer spectral index is produced by an AGN continuum 
which is getting brighter than the cold dust component of 
a spiral galaxy, we can interpret these plots saying again 
that the dominance of the AGN over the total power of 
a galaxy is well measured by its ionization indicators. 

Using our statistical tests, we have found that the 
[NeV] /[Nell] ratio correlates with the spectral index for 
all Seyfert galaxies (see the caption of Fig. 10 for the 
numerical results), while the [NcIII]/[NeII] ratio shows 
an increasing trend with the spectral index but not a 
statistical correlation. 

One advantage of the photometric indicators of AGN 
strength is that since they do not require spectroscopy, 
they will be much easier to measure for faint galaxies at 
high redshifts. 

5.2.3. Neon line ratios versus the [SIIIJ to [Sill] ratio 

In their study of the SINGS Spitzer legacy 
program. Dale et al. (2006) found that the 



[SIII]33.5/im/[SiII]34.8/im line ratio can be used to 
distinguish AGNs from starbursts. We therefore com- 
pared this ratio against the [NeV 14.3/im/[NeII]12.8/xm 
line ratio and against the [NeIII]15.5/um/ Nell]l2.8yum 
line ratio, to look for differences between the two Seyfert 
populations. We plotted in Fig. 11 (a, b) the two line 
ratio diagrams. Although of low statistical significance, 
it appears that the [SIII] to [Sill] ratio is somehow sepa- 
rating Seyfert I's from Seyfert 2's, in the sense that most 
of the former (8/12) have [SIII]33.5Aim/[SiII]34.8;um < 

0. 6, while most of the Seyfert 2's (8/10) have this ratio 
> 0.6. This threshold is also the one shown in the Dale 
et al. (2006) published diagram. However, the lack of 
correlation in Fig. 6 shows that the [SIII] /[Sill] ratio 
is less effective at picking out the strength of the AGN 
emission in detail. This is because it only spans a range 
of less than one order of magnitude, while the Ne line 
ratios span nearly two orders of magnitude. 

6. COMPARISON OF SEYFERT I'S AND 2'S 

The unification hypothesis states that Seyfert 1 and 
Seyfert 2 galaxies are intrinsically the same, but appear 
different at some wavelengths because of the much larger 
extinction along our lines-of-sight to the active nuclei 
in Seyfert 2's. At some infrared wavelength, the AGN 
emission becomes more isotropic, so that Seyfert 1 and 
Seyfert 2 of the same intrinsic luminosity should appear 
identical. Thus a fair test of the unification hypothe- 
sis requires matching Seyfert I's and 2's of comparable 
bolometric luminosity, and this is one reason we are an- 
alyzing the 12/im sample, which is relatively unbiased, 

1. e. equally sensitive to Seyfert Is or Seyfert 2s having 
the same bolometric fluxes. 

The equivalent widths of the PAH 11.25/xm or 
[NeII]12.8/xm emission features — which are unrelated 
to the AGN — are normalized by the underlying nuclear 
continuum at 12/im, which is predominantly coming from 
the active nucleus. This is why these EWs are anticor- 
related with the relative dominance of the AGN compo- 
nent. As expected, the EW of high-ionization lines such 
as [NeV] or [Nelll] is not correlated with AGN dominance 
because both those emission lines, and the underlying nu- 
clear continuum, is dominated by the AGN. The AGN 
itself produces [NeV] and [Nelll] that are about 0.03% of 
its bolometric power; since the 12/im continuum is 10% 
of the bolometric power for Seyfert galaxies this results 
in the "pure AGN" EWs of these fines of 0.3 /im. 

By normalizing the 12/xm continuum to other observ- 
ables which we believe are emitted roughly isotropically 
from AGN, we can compare directly the mid-IR contin- 
uum emitted from Seyfert I's and Seyfert 2's. Although 
it has long been known that Seyfert 2 nuclei arc relatively 
weaker in their 12/xm continuum (e.g. Edelson & Malkan 
1986; Edelson et al. 1987) the difference with Seyfert 
I's is not large in our sample. 

We observe that the average equivalent width (EQW) 
of the 11.25/im PAH for Seyfert I's is 0.38±0.39 /im while 
for Seyfert 2's it is 0.44±0.39 /J,m. These two values differ 
by less than la. This can be understood if some obscu- 
ration is still present in the mid-IR in Seyfert 2's. One 
popular possibility is that it is higher since the plane 
of the AGN torus is along our line of sight. But the 
effect must be quite weak. The EW's of the PAH at 
11.25/tm and that of the [NeII]12.8/tm are about 20% 
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larger in the Seyfert 2's, while the AGN-dominated [Ne 
III] and [NeV] EWs are also about 20% larger in Seyfert 
I's. Assuming the basic tenet of unification-that Seyfert 
I's and Seyfert 2's are intrinsically the same, this implies 
that only 20% of the 12/xm continuum from Seyfert 2's 
has been absorbed or blocked by geometric effects. This 
relatively modest implied optical depth to the 12/im con- 
tinuum of Seyfert 2's is far lower than the values required 
for "standard" thick dusty torus models (Pier & Krolik 
1992, 1993). However, recent models have allowed for 
the more likely possibility that the dust distribution is 
clumpy, so that some sight-lines to Seyfert 2's could in- 
tercept relatively little absorption (Elitzur 2006). This, 
combined with the likelihood of mid-IR continuum emis- 
sion from optically thin dust above the torus in both Sy 
I's and 2's, means that our data may not rule out these 
more sophisticated unification models. 

6.1. Host Galaxy ISM 

Photometric mid-IR studies of Edelson et al. (1987) 
Maiolino et al. (1998) indicated that Seyfert 2 galaxies 
more often have enhanced star formation than Seyfert 
I's. This finding has also been supported by the low res- 
olution IRS spectroscopy of Buchanan et al. (2006), 
who conclude that Seyfert 2 galaxies typically show 
stronger starburst contributions in their IR spectra than 
the Seyfert 1 galaxies. This result, if confirmed in a sam- 
ple which is not compromised by selection effects, would 
be in disagreement with the predictions of the unified 
scheme (Antonucci 1993). 

The PAH and H2 emission fluxes and luminosities, ac- 
cording to Rigopoulou et al. (2002), can discriminate be- 
tween starburst and AGN. In Fig. 12, we show the dia- 
gram with fluxes and luminosities (a and 6, respectively). 
The emission from small dust grains and warm molecules 
appears to be only correlated when the luminosities are 
taken, and not using the fluxes. Considering separately 
the Seyfert I's and 2's, we find similar correlations, there- 
fore there is no contradiction with the unification idea. 

We notice, however, that the correlation between lumi- 
nosities, that does not appear between the fluxes, could 
be induced by the distance factor and therefore it may 
not be real. 

6.2. Line ratios as density indicators 

Dudik et al. (2007) performed the first sys- 
tematic study of the [NeV](14.3/im/24.3/im) and 
[SIII](18.7/im/33.5/im) mid-infrared fine structure dou- 
blets of a heterogeneous sample of Seyferts (type I's 
and 2's), LINERs and quasars, from the Spitzer archive. 
They found that the [NeV] line ratio is often below the 
theoretical low-density limit, using a model that ignored 
absorption and stimulated emission. This appeared par- 
ticularly true for the type 2 AGNs. They claim that the 
ratios below the low-density limit are due to differential 
infrared extinction to the [NeV] emitting region caused 
by the obscuring torus, and that this ratio can trace the 
inclination angle of the torus to our line of sight. Be- 
cause the amount of extinction, the temperature of the 
emitting gas and the effects of absorption and stimulated 
emission are all unknown, they conclude that the [NeV] 
ratio cannot be used to infer the gas density. Similarly, 
they caution that the [SHI] line ratio is also useless be- 
cause mainly originated from the star forming regions of 



the host galaxies and the different slit aperture sizes used 
to observe the two lines prevent meaningful conclusions 
about the gas density in nearby galaxies. 

Even taking into account the caveats reported by 
Dudik et al. (2007), we nevertheless plotted in Fig. 13 
these two line ratios, which are in principle sensitive to 
the gas electron densities. The dashed lines show the low 
density limits for both ratios and the solid lines indicate 
the values of the density, as derived from the models of 
Dudik et al. (2007) assuming a temperature of 10,000K. 
Most of the sample galaxies lie close to the low density 
limits. This supports a thin spatially extended "Coronal 
Line Region" , rather than a very small " Intermediate 
Line Region" (Spinoglio & Malkan 1992). The densi- 
ties derived from the [NeV] ratio are of the order of 10^ 
cm~^. It appears that type 2 Seyfert's have /iZ(7/ier densi- 
ties of the highly ionized gas compared to type I's. Half 
of the Seyfert 2's have Ue > lO'^'^ cm^'^, while all Seyfert 
I's have densities lower than this value. The densities 
derived from the [SHI] ratio are slightly lower, clustering 
around Ug ~ 10^ cm~^. 

We do not find the same results as reported from Dudik 
et al. (2007): only two Seyfert 2 galaxies (NGC1125 
and NGC3660) are below the low density limit derived 
for the [NeV] line ratio, and the same number of Seyfert 
I's (ESO012-G021 and NGC4748, that becomes four in- 
cluding the objects having upper limits, namely MRK817 
and NGC4602). We notice, however, a greater number of 
galaxies falling below the [SHI] line ratio density limit. 
While the [NeV] line ratio cannot be affected by aper- 
ture effects, because the [NeV] emission surely originates 
from the AGN, we cannot rule out the presence of these 
effects for the [SHI] line ratio, which is excited also in 
HII regions. About this latter, however, we can exclude 
such effects in 11 galaxies (about 1/3 of the objects), 
for which we have the [SIII]18.71/Ltm line measured in 
the LH slit (sec Table 2). For the other objects with a 
low [SIII](18.7/um/33.5/um) ratio, we have used the ex- 
tendness factors (Table 4) to correct the [SIII] 17.81 jim 
line flux, assuming that increasing the aperture will in- 
crease by the same proportion continuum and line fluxes. 
Because the extendness factors were defined by measur- 
ing the continuum around the wavelength of 19 /im, we 
basically assume that the equivalent width of the [SIII] 
line does not change with the aperture. We show in the 
diagram of Fig. 13, both the measured and the aperture- 
corrected [SIII] line ratios directly above them. By ap- 
plying this simple correction, we find that only 4 galaxies 
(IRAS00198-7926 and NGC513, which have' been mea- 
sured with the LH slit, and ESO141-G055 and NGC6890) 
have a [SIII] ratio below low density limit. 

7. CONCLUSIONS 

We have identified a family of IRS obscrvables which 
quantify the proportion of the total IR emission coming 
from a Seyfert nucleus, all of which are intercorrelated 
with each other. The ratios of ionic fine structure lines 
[NeV] /[Nell] and [OIV]/[NeII] were aheady proposed to 
measure the importance of the AGN component. We also 
see that [NeIII]/[NeII] and OIV or [NeV]/[SiII]35/xm can 
be used to quantify the AGN dominance. 

It was also known from ISO that the equivalent width 
of the 6.2/xm (and 7.7/L(m) PAH features is inversely re- 
lated to the AGN dominance; we find that the same holds 



Spitzer spectra of Seyfert galaxies 



for the equivalent width of the 11.25/im PAH feature. 
We also discovered two additional IRS observablcs: the 
equivalent width of [NeII]12.8/Ltm and the extendness of 
the IQyum continuum, which also quantify the dominance 
of the AGN component compared with the emission from 
the underlying spiral galaxy. All of these observables are 
correlated with each other, since they are measuring this 
same astrophysical quantity, and they are all correlated 
with the hardness of the far-IR continuum, since a more 
dominant AGN component is already known to be cor- 
related with hotter dust. 

There is no clear indication that recent star formation 
is much more important on average in the Seyfert 2's in 
our sample, compared with the Seyfert I's. Although 
the Seyfert I's generally tend to have more dominant 
AGN than the Seyfert 2's, there is a strong overlap be- 
tween these two classes. The relatively small difference 
between the averages of these observables for Seyfert I's 
and Seyfert 2's indicates that those two AGN categories 
are not extremely different in the mid-IR range. Thus for 
example, the AGN in Seyfert 2's are clearly observable 
in their 10/im contimmm emission, and in [NcIII], [OIV] 
and [NeV], to almost the same degree as in Seyfert I's, 
which is not consistent with some proposed torus models 
for Seyfert 1/2 unification. 
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Finally, we do not fully confirm the observational 
claims of Dudik et al. (2007) : 75 - 80% of our Seyfert 
I's and Seyfert 2's show [SHI] or [NeV] densities larger 
than the low-density limit. In fact, after applying plausi- 
ble aperture corrections to the [SIII] line ratio, only three 
Seyfert 2's and one Seyfert 1 have a [SIII] line ratio be- 
low the density limit. A similar result is also found for 
the [NeV] ratio, for which we have two Seyfert I's (that 
became four including the upper limits) and 2 Seyfert 2's 
below the low density limit. These few cases do not in 
our view require enormous dust extinction values. 
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APPENDIX 

ADDITIONAL OBSERVATIONAL MEASUREMENTS 

Weak ionic fine structure lines 

In the IRS high resolution spectra presented, besides the common lines reported in Table 2 and 3, we have also mea- 
sured a few fine structure lines in the spectra of a few galaxies. These include the [Aril] 13. 10/im and the [MgV]13.52/im 
lines, as well as five [Fell] lines and the [FeIII]22.93/im line. Among the [Fell] lines, we notice that the line at 25.99/xm 
can be blended with the much brighter 25.89/;im [OIV] line, because the spectral resolution is only about half of the 
lines separation. Except in the few cases where the lines were well separated, we did not attempt to deblend them, 
but we only gave a measurement of the [OIV] line. We refer to a future paper, where we will compare the data of our 
galaxies with photoionization models, for more accurate measurements of the [Fell] 25.99/im line. We report in table 
Al the fluxes, derived from gaussian fits, of the weak lines detected. 

Equivalent widths of the key fine structure lines 

We report in Table A2 the equivalent widths (in /xm) of the brightest lines, which have been used in the diagrams 
and correlations presented in the main text. 
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TABLE Al 
Journal of Spitzer IRS observations 





R.A. 


Dec. 












SH 


LH 




NAME 


(J2000.0) 


(J2000.0) 


type 


z 


Fl2um 


F25um 


Obs. Date 


Int. time 


Int. time 






h m s 


deg "' 






(Jy) 


(Jy) 




(sec.) 


(sec.) 




IRAS00198-7926 


00:21:57.0 


-79:10:14 


Sy2 


0.072800 


0.28 


1.15 


2006/06/28 


4x30 


2x60 




ESO012-G021 


00:40:47.8 


-79:14:27 


Syl 


0.032840 


0.17 


0.25 


2006/06/28 


4x30 


2x60 




IRAS00521-7054 


00:53:56.2 


-70:38:03 


Sy2 


0.068900 


0.28 


0.80 


2006/06/29 


4x30 


2x60 




ESO541-IG012 


01:02:17.5 


-19:40:09 


Sy2 


0.056552 


0.20 


0.36 


2006/08/03 


4x30 


2x60 


bzer 


NGC0424 


01:11:27.5 


-38:05:01 


Sy2 


0.011764 


1.10 


1.73 


2006/06/28 


2x30 


4x 14 


NGC0513 


01:24:26.8 


-1-33:47:58 


Sy2 


0.019544 


0.17 


0.28 


2006/09/17 


2x30 


4x60 


spe( 


MRK1034 NEDOl 


02:23:18.9 


-1-32:11:18 


Syl 


0.033633 


0.25 


0.69 


2006/09/19 


2x30 


4x14 


MRK1034 NED02 


02:23:22.0 


4-32:11:49 


Syl 


0.033710 






2006/09/19 


2x30 


4x14 


y- 


ESO545-G013 


02:24:40.2 


-19:08:27 


Syl 


0.033730 


0.17 


0.35 


2006/09/08 


4x30 


2x60 


i-i 








oyi 


U.UIDDOZ 


U.DJ. 


1 

J..01 


onnf^ /HQ /I Q 


Z A OU 




O 


NGC1125 


02:51:40.4 


-16:39:02 


Sy2 


0.010931 


0.17 


0.83 


2006/09/08 


4x30 


2x60 




ESO033-G002 


04:55:59.6 


-75:32:27 


Sy2 


0.018100 


0.20 


0.45 


2006/07/27 


4x30 


2x60 


MRK0006 


06:52:12.2 


-1-74:25:37 


Syl 


0.018813 


0.22 


0.69 


2006/11/16 


4x30 


2x60 


(D 
i-i 
ci- 


MRK0009 


07:36:57.0 


-1-58:46:13 


Syl 


0.039874 


0.21 


0.44 


2006/11/16 


4x30 


2x60 


NGC3516 


11:06:47.5 


-1-72:34:07 


Svl 


0.008836 


0.43 


0.89 


2006/11/16 


2x30 


4x14 




NGC3660 


11:23:32.2 


-08:39:30 


Sy2 


0.012285 


0.19 


0.22 


2006/06/27 


2x30 


4x14 


P 


NGC4501 


12:31:59.0 


-1-14:25:10 


Sy2 


0.007609 


1.02 


1.28 


2006/06/27 


2x30 


4x14 


? 


NGC4602 


12:40:36.5 


-05:07:55 


Syl 


0.008469 


0.54 


0.57 


2006/06/27 


2x30 


4x14 


CD 


TOLOL01238-364 


12:40:52.9 


-36:45:22 


Sy2 


0.010924 


0.64 


2.26 


2006/07/03 


2x30 


4x14 




NGC4748 


12:52:12.4 


-13:24:54 


Syl 


0.014630 


0.17 


0.37 


2006/07/01 


4x30 


2x60 




NGC4968 


13:07:06.0 


-23:40:43 


Sy2 


0.009863 


0.39 


1.05 


2006/06/29 


2x30 


4x14 




MCG-03-34-063 


13:22:19.0 


-16:42:30 




0.021328 


0.95 


2.88 


2006/06/29 


2x30 


4x14 




MRK0817 


14:36:22.1 


58:47:39 


Syl 


0.031455 


0.33 


1.18 


2007/03/20 


4x30 


2x60 




IRASF15091-2107 


15:11:59.8 


-21:19:02 


Syl 


0.044607 


0.23 


0.50 


2007/03/08 


4x30 


2x60 




ESO141-G055 


19:21:14.1 


-58:40:13 


Syl 


0.036000 


0.24 


0.35 


2007/04/30 


4x30 


2x60 




NGC6890 


20:18:18.1 


-44:48:23 


Sy2 


0.008069 


0.34 


0.65 


2006/10/25 


4x30 


2x60 




MRK0897 


21:07:45.8 


-1-03:52:40 


Sy2 


0.026340 


0.37 


0.86 


2006/11/16 


4x30 


2x60 




IRASF22017+0319 


22:04:19.2 


-1-03:33:50 


Sy2 


0.061100 


0.29 


0.72 


2006/11/19 


4x30 


2x60 




NGC7496 


23:09:47.2 


-43:25:40 




0.005500 


0.35 


1.60 


2006/06/28 


2x30 


4x14 




NGC7590 


23:18:55.0 


-42:14:17 




0.005255 


0.56 


0.83 


2006/06/28 


2x30 


4x14 





CD 



TABLE A2 
Common ionic fine structure lines 



Line fluxes (10 erg s ^cm ^) in SH Line fluxes (10 erg s ^cm ^) in LH 
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[SIV] 


[Nell] 


[NcV] 
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o oi -Un OA 


Z.UZ±U.o4 




1 fic:_i_n on 
l.DoitU.zU 


i. iOniU.oo 


A Ofi-Un 

4,yozizU.O 




A fif^-un 0/1 
4.DDltU.y4 




Q OQ-Un A A 

o.yoiLU.44 


Q '7n_i_n QQ 
o. i U±U.v>o 


1 1 n_Ln 7 
iD.lUitU. / 


io.40±U.4z 


D.yDltU.Oi 




0.v5 ( itU.4o 


OK Q -1- n K'\ 

zo.odiU.oi 


O QO-UI OQ 

Vi.oZzti.zo 


Q 1 /I _L1 /II 
o. i4iti.4i 


IN OOUOio 


O 'WJ-Oi AO 
Z. 1 1 ztU.4Z 


1 o 7fi-J-n oi 


1 OI -un Q/1 


4.4oztU.Zy 


D. i DitU.OZ 


D.Z I itu. io 


1 no-un no 
i.uyitu.uy 


fi. KA-\-n /10 
D.04ztU.4y 


i4.0U±U.O 1 


ov /in-un c^fi 
z ( .4y±:U.OD 


lVir\,l\-lUo4 INJl/JJUl 




lo.OoztU.oz 


< i.i 


1.4 ( ±U.oU 


Q i^OJ_n A Q 




^ n c;n 

< u.oy 


^ n '71 

< U. M 


"in onj-n qi 
iU.zUzbU.oi 


0/1 A 1 _uo nn 
Z4.4iitz.UU 


iVirVxvlUo^ iNrLjUUZ 


1 QQ_|_n on 


OA QQ_|_n Of: 


1 .UDlCU.oO 




Q nv-i-fi 
y .u / mu.o^ 




U.DU 


Z.DOICU.O ( 


1 n -1-1 1 zl 


r:i 71 j_o nn 


ESO545-G013 


5.32±0.22 


10.05±0.16 


4.80±0.25 


10.50±0.30 


3.59±0.41 


5.87±0.13 


3.22±0.11 


11.54±0.15 


8.31±0.45 


20.00±0.54 


NGC0931 


10.70±0.48 


5.47±0.43 


14.30±0.78 


15.41±0.43 


4.86±0.61 




13.67±0.43 


42.60±0.38 


11.97±1.28 


13.72±1.70 


NGC1125 


6.07±0.22 


16.37±0.20 


5.09±0.26 


15.55±0.23 


10.99±0.27 




9.69±0.31 


40.36±0.53 


23.94±0.97 


31.32±0.78 


ESO033-G002 


5.54±0.27 


2.13±0.19 


6.34±0.16 


9.22±0.18 


3.99±0.29 




5.27±0.17 


14.39±0.18 


2.79±0.65 


4.53±0.44 


MRK0006 


16.69±0.36 


28.00±0.23 


9.39±0.19 


49.34±0.32 


14.10±0.25 




10.43±0.21 


48.24±0.27 


14.09±0.54 


36.40±0.71 


MRK0009 


2.37±0.23 


3.23±0.20 


2.21±0.22 


1.90±0.20 


2.38±0.50 


3.44±0.15 


2.24 ±0.19 


5.55±0.21 


3.94±0.49 


7.32±0.77 


NGC3516 


13.33±0.38 


8.07±0.25 


7.88±0.50 


17.72±0.33 


5.86±0.35 




10.39±0.33 


46.92±0.35 


9.52±0.96 


22.14±0.54 


NGC3660 


1.48±0.22 


6.51±0.23 


0.98±0.23 


1.49±0.18 


3.66±0.28 




1.66±0.20 


3.61±0.25 


9.52±1.14 


9.54±1.10 


NGC4501 


< 3. 


7.02±0.27 


< 1.5 


4.72±0.29 


1.97±0.31 




< 3.6 


4.22±0.28 


7.9 ±0.6 


16.70±0.87 


NGC4602 


< 1.2 


7.57±0.25 


0.82±0.09 


0.63±0.19 


3.20±0.23 




< 1.2 


< 2.3 


8.68±0.48 


15.30±0.83 


TOLOL01238-364 


5.70±0.28 


45.15±0.54 


11.15±0.61 


27.00±0.44 


16.32±0.48 




5.35±0.67 


21.21±0.58 


32.80±1.62 


44.99±1.62 


NGC4748 


9.87±0.18 


7.37±0.19 


6.75±0.19 


15.93±0.20 


7.61±0.22 


16.33±0.47 


20.00 ±0.17 


81.93±0.21 


20.42±0.68 


29.46±0.61 


NGC4968 


9.63±0.37 


24.90±0.30 


17.57±0.57 


33.80±0.44 


15.10±0.58 




10.60±0.55 


33.70±0.45 


14.90±1.14 


22.36±1.15 


MCG-03-34-063 


< 1.2 


6.02±0.27 


< 2.4 


< 3. 


3.11±0.36 




< 3. 


< 3. 


6.67±0.65 


6.71±1.60 


MRK0817 


1.53±0.38 


3.83±0.22 


1.86±0.19 


4.58±0.41 


2.76±0.60 




3.60±0.44 


6.53±0.29 


< 3.21 


< 3.45 


IRASF15091-2107 


7.29±0.22 


11.52±0.15 


8.48±0.30 


16.29±0.29 




9.95±0.15 


8.12±0.15 


31.03±0.30 


12.26±0.62 


12.27±1.49 


ESO141-G055 


3.45±0.22 


2.24±0.21 


2.25±0.24 


5.62±0.27 


1.75±0.52 




1.62±0.12 


7.26±0.12 


5.45±0.50 


8.85±0.82 


NGC6890 


2.92±0.28 


11.32±0.32 


5.77±0.29 


6.57 ±0.14 


4.34±0.29 




3.77±0.18 


10.10 ±0.17 


16.97±0.72 


26.54±0.77 


MRK0897 


< 1.5 


24.03±0.19 


1.06±0.21 


4.38 ±0.13 


14.91±0.29 


15.52±0.24 


< 0.8 


0.62±0.19 


22.30±0.70 


21.44±0.85 


IRASF22017+0319 


10.31±0.25 


5.95±0.37 


8.33±0.26 


14.07 ±0.29 




6.10 ±0.25 


9.40±0.30 


29.04±0.30 


9.33±0.44 




NGC7496 


1.3 ±0.4 


48.08±0.37 


< 1.8 


6.67±0.35 


23.48±0.50 




< 2.4 


< 2.4 


39.47±1.17 


44.58±2.11 


NGC7590 


< 3 


7.78±0.20 


< 1.5 


3.49±0.27 


5.37±0.37 




< 1.2 


5.60±0.25 


15.69±0.83 


26.50±0.95 
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TABLE A3 

Molecular Hydrogen lines and PAH emission feature at 11.25/um 



Line fluxes (10 erg s ^crn 





rl2 ^\-~>] 


112 ^\^} 


112 '^{^j 


112 ^{^) 


PAH 


r^v.^. vv . 






llZ.ZO/^//t 1 






(-\-\ OKiim^ 


(pm) 








(•^) 






(a) 


IRAS00198-7926 


2.31±0.51 


4.22±0.36 


3.92±0.29 


< 1.47 


24.7 


-0.117 


ESO012-G021 


2.35±0.14 


1.47±0.13 


3.93±0.34 


< 0.54 


65.5 


-0.404 


IRAS00521-7054 


2.75±0.30 


1.78±0.30 


2.40±0.21 


< 1.11 






ESO541-IG012 


< 0.65 


< 0.72 


1.24±0.19 


< 1.2 


86.4 


-0.434 


NGC0424 




< 1.14 


2.01±0.40 


< 1.53 


11.0 


-0.021 


NGC0513 




< 8.60 


2.67±0.29 


1.54±0.10 


63.0 


-0.704 


MRK1034 NEDOl 


< 1..5 


< 0.96 


3.50±0.36 


< 2.4 


68.8 


-1.054 


MRK1034 NED02 


3.14±0.32 


< 0.72 


6.79±0.56 


< 1.11 


119.0 


-0.754 


ESO545-G013 


1.03±0.20 


< 0.48 


1.57±0.36 


1.01±0.18 


20.3 


-0.494 


NGC0931 




< 1.29 


< 1.83 


< 1.14 


24.6 


-0.064 


NGC1125 




< 0.60 


1.28±0.32 


< 1.53 


19.9 


-0.655 


ESO033-G002 




1.29±0.22 


2.44±0.25 


0.82±0.19 


12.8 


-0.080 


MRK0006 




1.10±0.24 


4.26±0.26 


1.49±0.29 


18.1 


-0.092 


MRK0009 


1.33±0.27 


< 0.53 


1.19±0.32 


< 0.63 


13.1 


-0.078 


NGC3516 




1.00±0.39 


3.93±0.41 


< 1.05 


22.1 


-0.076 


NGC3660 




< 0.69 


2.13±0.24 


< 0.75 


21.7 


-0.635 


NGC4501 




4.35±0.28 


10.07±0.29 


2.71±0.37 


30.2 


-0.997 


NGC4602 




1.84±0.19 


2.44±0.23 


2.94±0.26 


24.1 


-1.190 


TOLOL01238-364 




< 1.62 


a 44-1-0 46 


< 1.74 


94.3 


-0.201 


NGC4748 




< 0.57 


2.21±0.27 


0.71±0.22 


18.6 


-0.615 


NGC4968 




< 1 


3.34±0.39 


< 1.35 


69.3 


-0.203 


MCG -03-34-063 




< 0.5 


2.04±0.30 


< 1.33 


10.8 


-1.043 


MRK0817 




< 0.6 


1.55±0.45 


< 1.46 


15.8 


-0.242 


IRASF15091-2107 


2.06±0.24 


1.40±0.21 


4.04±0.24 


< 0.84 


28.4 


-0.161 


ESO141-055 




< 0.6 


< 1.09 


< 0.74 


16.2 


-0.084 


NGC6890 




< 0.96 


1.31±0.32 


1.68±0.20 


40.6 


-0.382 


MRK0897 


< 0.57 


1.07±0.09 


2.84±0.24 


0.98±0.18 


78.7 


-1.279 


IRASF22017-t-0319 


1.84±0.24 


0.88±0.25 


1.90±0.27 


< 0.76 


24.0 


-0.072 


NGC7496 




< 1.11 


4.95±0.33 


< 1.70 


105.0 


-0.161 


NGC7590 




1.63±0.25 


2.64±0.28 


2.36±0.33 


38.8 


-0.081 


* Columns (1) to (4) 


give the line fluxes of the 


H2 rotational 


lines in units 


of 10""^* erg 





Column (5) and (6) give the flux and the equivalent width of the PAH 11.25 iJ.m emission feature, 
respectively. 
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TABLE A4 

Photometric fluxes and extendness factors 



NAME 


Fi9Mm(SH) 


Fl9.5um(LH) 

(Jy) 


R 


Class 




(Jy) 






TRA'^nniQS 7Q9fi 


n c;74_|_n 099 


79Q-I-0 01 ^ 


1.27 


jj 




U. -L'±U_L.U.U -LU 


1 744-0 OOQ 


1.19 


jj 


TRA^nn'^91 7V\^A 


'^074-0 091 


'^'^1 4-0 01 


1.05 


J 


J_lJ O v_7 rt -L -L v_T V X Zi 


A ocn-l-n ni 9 


'iofi-i-o 00s 

u.ouu^u. uuo 


1.09 


J 


i\rpr^nzi9zi 


1 on 1 _|_n 09*^ 


1 CI 7J^o 09*^ 

L.OL 1 3ZU.UZO 


1 OQ 


J 




n 1 09-1-0 091 


1 f^O-l-0 00'^ 


1.47 


JJ 


IVll ViV -LU04l in llJLJy} X. 


nfifi-l-n 01 ^ 


0SQ4-0 007 


1.34 


JJ 


IVll V,iV _LU^ J i- i-N HjL^kjZj 


9904-0 090 


Qi 04-0 01 1 


1.41 


JJ 




1 07-1-0 01 


1 404-0 OOfi 


1.31 


JJ 




(\(\AA-(\ 094 


QqQ4-o 01 1 


1 26 


JJ 


NOfl 1 9^ 


A 9'3c;_|_o 01 ZL 


44^^4-0 01 1 


1.91 


III 




A QCC_|_n Al 9 


A QS74-n 007 


1.09 


J 


IVl X\, Xv u u u u 




A c;77_|_n Al 9 


1.00 


J 


MRKnnnQ 

IVX X\, JVU U U i7 


A Q99_|_n A14 


A Qc;n4-n 007 


1.09 


J 


IN '^^OO -LU 




A K79-|-n 01 Q 


1 14 


J 




A nfi7-i-n m 9 


CiR7-\-C\ 01 9 


1.30 


JJ 


il V_T V_y '-t 1 


n n'^2-i-n ni 

\j •yjtj^-L.yj ^yj 


A nc4_|_n Al 9 


1.69 


III 




049-1-0 OOS 


1 01 4-0 01 9 


2.40 


III 


T0T,0T,01 2SS-Sfi4 


1 9Q0-I-O 04*^ 


1 'i22+n 021 

i- . '.J Z^ 1 . U Z i 


1.18 


JJ 


NGC4748 


0.089±0.010 


0.223±0.008 


2.50 


III 


NGC4968 


0.812±0.029 


0.817±0.014 


1.00 


I 


MCG-03-34-063 


0.020±0.010 


0.042±0.009 


2.09 


III 


MRK817 


0.670±0.016 


0.778±0.014 


1.16 


II 


IRASF15091-2107 


0.337±0.012 


0.388±0.008 


1.15 


I 


ESO141-055 


0.249±0.013 


0.288±0.008 


1.16 


II 


NGC6890 


0.328±0.019 


0.381±0.007 


1.16 


II 


MRK0897 


0.120±0.010 


0.164±0.006 


1.37 


II 


IRASF22017+0319 


0.533±0.020 


0.583±0.010 


1.09 


I 


NGC7496 


0.658±0.031 


0.806±0.013 


1.22 


II 


NGC7590 


0.032±0.010 


0.096±0.011 


3.00 


III 



TABLE A5 

H2 Estimated Temperatures and Masses 



NAME 


T(O-l) 


T(l-2) 


T(2-3) 


H2 Mass 




(K) 


(K) 


(K) 


(106 Mq) 


IRAS00198-7926 


> 229 


573 


241 


< 42.3 


ESO012-G021 


> 422 


300 


351 


< 4.19 


IRAS00521-7054 


> 209 


434 


346 


< 26.7 


ESO541-IG012 


> 160 


< 374 




< 4.37 


NGC0424 


> 173 


< 370 




< 0.9 


NGC0513 


191 






0.88 


MRK1034 NEDOl 


182 


< 264 




11.8 


MRK1034 NED02 


> 173 


< 192 


> 185 


< 25.6 


ESO545-G013 


184 


< 275 


> 402 


5.18 


NGC1125 


> 150 


< 334 




< 0.73 


ESO033-G002 


241 


356 




1.4 


MRK0006 


2.36 


257 




< 2.8 


MRK0009 


> 198 


< 325 


> 435 


< 4.8 


NGC3516 


> 271 


256 




< 0.48 


NGC3660 


> 137 


< 282 




2.0 


NGC4501 


270 


321 




0.91 


NGC4602 


150 


439 




0.84 


TOLOL01238-364 


> 202 


< 335 




< 0.97 


NGC4748 


246 


< 257 




0.82 


NGC4968 


> 222 


< 273 




< 0.65 


MCG-03-34^063 


> 183 


< 252 




< 2.7 


MRK817 


> 162 


< 305 




< 0.28 


IRASF15091-2107 


> 198 


290 


339 


< 3.8 


NGC6890 


147 


< 431 




0.43 


MRK0897 


238 


301 


< 239 


3.6 


IRASF22017+0319 


> 223 


332 


396 


< 4.0 


NGC7496 


> 239 


< 244 




< 13.5 


NGC7590 


164 


388 




1.0 
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TABLE A6 

Table A1 - Rare ionic fine structure lines 






Line fluxes (10 ^'^ erg s ^cm 








NAME 


[ArV] 
(13.10//m) 


[MgV] [Fell] [Felll] [Fell] 
(13.52/nm) (17.94;im) (22.93;im) (24.52/im) 


[Fell] 
(25.99/im) 


[Fell] 
(35.35//m) 


[Fell] 
(35.77//m) 


NGC0513 

MRK1034 NEDOl 

MRK1034 NED02 

ESO545-G013 

NGC1125 

MRK0006 

MRK0009 

NGC3660 

NGC4501 




0.35±0.10 

2.39±0.28 

2.07±0.15 

1.40±0.20 


1.41±0.07 
1.04±0.27 
1.95±0.35 

2.62±0.28 


8.09±0.77 





TOLOL01238-364 ••• 2.59±0.24 

NGC4968 ■•• 2.44±0.22 

MCG-03-34-063 ■•• ••■ ••■ ••• ■•■ 9.10±1.48 

MRK0817 ••• ••■ ••• •■• < 1.0 ■•■ 14.10±1.25 
IRASF15091-2107 1.37±0.19 < 0.96 

MRK0897 ••■ ■■• ■•• 1.27±0.12 ••• 1.49±0.19 

NGr7196 ••• ••• ••• ••• ••• ••• ••• 5.72±1.8S 



TABLE A7 

Table A2 - Equivalent widths of key lines 



NAME [NeII]12.81/im [NeV]14.32^tm [NeIII]15.56/im H2 17.02^jm [OIV]25.89Aim 





(/im) 


(urn) 


(//m) 


(Mm) 


(/urn) 


IRAS00198-7926 


-0.009 


-0.023 


-0.029 


-0.009 


-0.065 


ESO012-G021 


-0.055 


-0.023 


-0.048 


-0.032 


-0.159 


IRAS00521-7054 


-0.008 


-0.010 


-0.016 


-0.006 


-0.027 


ESO541-IG012 


-0.006 


-0.008 


-0.008 


-0.004 


-0.043 


NGC0424 


-0.005 


-0.011 


-0.013 


-0.003 


-0.044 


NGC0513 


-0.072 


-0.017 


-0.049 


-0.023 


-0.072 


MRK1034 NEDOl 


-0.129 




-0.025 


-0.041 




MRK1034 NED02 


-0.107 


-0.012 


-0.018 


-0.025 


-0.014 


ESO545-G013 


-0.108 


-0.064 


-0.119 


-0.017 


-0.178 


NGC0931 


-0.007 


-0.017 


-0.024 




-0.112 


NGC1125 


-0.077 


-0.027 


-0.090 


-0.004 


-0.104 


ESO033-G002 


-0.003 


-0.018 


-0.028 


-0.008 


-0.098 


MRK0006 


-0.064 


-0.020 


-0.100 


-0.009 


-0.190 


MRK0009 


-0.010 


-0.007 


-0.006 


-0.004 


-0.035 


NGC3516 


-0.013 


-0.017 


-0.033 


-0.008 


-0.129 


NGC3660 


-0.095 


-0.015 


-0.029 


-0.033 


-0.065 


NGC4501 


-0.107 




-0.215 


-0.248 


-0.164 


NGC4602 


-0.122 


-0.027 


-0.043 


-0.056 




TOLOL01238-364 


-0.047 


-0.013 


-0.029 


-0.004 


-0.025 


NGC4748 


-0.096 


-0.113 


-0.248 


-0.034 


-0.658 


NGC4968 


-0.032 


-0.027 


-0.048 


-0.005 


-0.084 


MCG-03-34-063 


-0.162 






-0.126 




MRK0817 


-0.006 


-0.003 


-0.007 


-0.002 


-0.016 


IRASF15091-2107 


-0.036 


-0.034 


-0.057 


-0.013 


-0.159 


ESO141-G055 


-0.007 


-0.010 


-0.023 




-0.064 


NGC6890 


-0.038 


-0.021 


-0.024 


-0.005 


-0.043 


MRK0897 


-0.128 


-0.012 


-0.057 


-0.026 


-0.004 


IRASF22017-I-0319 


-0.015 


-0.020 


-0.037 


-0.005 


-0.121 


NGC7496 


-0.089 




-0.016 


-0.009 




NGC7590 


-0.098 




-0.113 


-0.061 


-0.090 
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IRAS00198-7926 



S(3) 


S(2) 


S(l) 






H3 



H.^ S(3) 9 



[Sni]l8.7^fm 





[Ntiilll 


[NeV] [SIlIl 1 


[jNeV] 


[SHI] 


[oiv] 



20 

\(rest) (,im) 



ES0012-G021 



1^ 



#11 



[NeV]H,3H/(m 



^1 



[NeV] [Sill] 



[OIV] 



[Sill] 



20 

A(rest) (,im) 



IRAS0052 1-7054 



S(3) 


S(2) 


s(i) 









Ha 3(2) 9.66fim [SIVl 1 0.5 l^j.m 



;i«^' 



/ HgS(2) 



|NeVl24.32/jm 

A 



JNbIII 


[Ntillll 


[NeV] ]Sin] 1 


[NeV] 


[Sill] 


[OIV] [ 



10 




20 






A(re3t.) (//m) 






NGC0424 




5(3) 




s(l) 






H3 1 



FS lines - 
"aiii]- 

30 



S(l) 17.04/^rr 



Si 
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Fig. Al. — Spitzer IRS SH and LH specra of the observed Seyfert galaxies. Wavelengths have been shifted to the galaxies rest frames. 
For the fainter lines (those with a S/N < 10.) the blow-ups with their theoretical wavelengths are also given. 
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Fig. A2. — H2 excitation diagrams. For each measured line, the natural logarithm of the level population normalized to its statistical 
weight is plotted against the upper level energy (in temperature units). For each pair of adjacent transitions the connecting line is shown, 
whose inverse value represents the gas temperature: the dashed line connects the S(0) and S(l) detections, the dotted line the S{1) and 
S(2), the solid line the S(2) and S{3). Upper limits have been used to obtain limiting slopes and hence limiting temperatures and masses 
(see text). 
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Fig. A3. — a:[NeV]14.3/xm/[NeII]12.8/im line ratio versus the PAH 11.25/im equivalent width. Solid triangles show the good detections 
of Seyfert type 2 galaxies and solid squares of Seyfert I's, open triangles and open squares show upper limits of Seyfert 2's and I's, 
respectively. The circles indicate MCG-3-34-63, NGC7496 and NGC7590, which are not classified as a Seyfert galaxies (see text), b: 
[OIV]14.3/^m/[NeII]12.8/xm line ratio versus the PAH 11.25/jm equivalent width. Symbols and lines as in the previous figure. 




-1.5 



-1 -0.5 

EW PAII (11.25Mm) 




iJW PAH (il.£5;/m) 



Fig. A4. — a:[NeV]14.3^m/[SiII]34.8/im line ratio versus the PAH 11.25/im equivalent width. Symbols and lines as in the previous figure, 
b: [NeHI] 15.5^m/[NeH] 12.8/im line ratio versus the PAH 11.25/im equivalent width. Symbols and lines as in the previous figure. 
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Fig. A5. — [NeII]12.8/jm equivalent width versus the PAH equivalent width. Symbols in the previous figure. We note that, for graphical 
reasons, in this diagram and in the following ones in which the equivalent widths, covering a large range, are shown, the logarithm of the 
inverse of the actual equivalent width is plotted. 
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Fig. A6. — PAH11.25/im and [NeII]12.8/im equivalent widths versus extendness parameter. Extendness of 1.0 corresponds to completely 
unresolved (point-like) 19/^m continuum, which is therefore completely AGN-dominated. In both figures the non-linear effects of extended 
emission from the host galaxy become important for extendness parameters above 1.3. Symbols as in the previous figure. 
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Fig. A7. — [NeIII]15.5/xm line equivalent width versus extendness parameter. Symbols as in the previous figure. 
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Fig. A8. — [NeV]14.3/im and [OIV]25.9^m line equivalent widths versus extendness parameter. As expected, there are no significant 
correlations, because the equivalent width is the ratio of high-ionization emission to the hot dust continuum, and both of these arise mostly 
from the AGN. Symbols as in the previous figure. 
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Fig. A9. — [NeV]14.3/im/[NeII]12.8/^m line ratio versus 12/im luminosity. 
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Fig. AlO. — a: [NeV]14.3/im/[NeII]12.8/^m line ratio versus the 25-60 ^tm spectral index. A correlation is present for all Seyfert's data, 
with a slope of 1.33 using the weighted least squares fit and 1.32 ± 0.77 (1 cr) using the bootstrap method, b: [NeIII]15.5/im/[NeII]12.8/im 
line ratio versus the 25-60 /im spectral index. 
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Fig. All. — a:[NeV]14.3/.tm/[NeII]12.8/im line ratio versus the [SIII]33.5/.tm/[SiII]34.8/.tm line ratio. Symbols as in the previous figures, 
b: [NeIII]15.5/im/[NeII]12.8/xm line ratio versus the [SIII]33.5/im/[SiII]34.8/im line ratio. Symbols as in the previous figures. Although 
uncorrelated with ionization level, the [SIII]/[SiII] ratio does roughly separate the Seyfert I's and 2's, although with substantial overlap. 




Fig. A12. — a: H2 17.04 line fiux versus PAH 11.25^m integrated flux. No correlation is present, b: H2 17.04 /im line luminosity 
versus PAH 11.25/.tm luminosity. The solid line is the least squares fit to all the Seyfert's data with a slope of 1.16, using a weighted least 
squares fit and 1.14 ± 0.13, using the bootstrap method. 
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Fig. A13. — [NeV]14.3/im/[NeV]24.3/im line ratio versus the [SIII]18.7^m/[SIII]33.5/im line ratio. Symbols as in the previous figures, 
except for the fact that we indicate here as open crossed symbols those for which the [SIII]18.71/im line was measured only with the small 
(SH) aperture and the above filled symbols directly above them show the corrected [SIII] line ratio for that object. The dashed lines show 
the low density limits (see text). The solid lines at the top and at the right give the corresponding electron densities. 



